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"Polaroid  sets  and  functions  have  been  introduced  as  a  new  tool,  with 
applications  in  non-linear  programming,  particularly  in  quasi-concave  and 
Integer  optimization  problems  over  a  linearly  constrained  set  of  feasible 
solutions . 

The  name  polar  programming  applies  to  a  general  class  of  non-linear 
mathematical  programming  problems  which  can  be  solved  by  the  polaroid 
approach.  In  integer  programming  polaroids  yield  non-trivial  extensions  of 
the  intersection  cut  approach. 


!fh4 «  paper  builds  on  the  properties  of  polaroid  sets  (particularly  com¬ 
plete  convex  polaroids)  and  focuses  on  the  following  intersection  problem: 

Given  a  point  x  belonging  to  the  polaroid  set  P*,  find  the  inter¬ 
section  point  u*  of  a  one-dimensional  ray  u  with  the  boundary  of  P*  a 
u*  c  (bd  F*  flu). 

We  also  present  a  theoretical  comparison  of  the  relative  merits  of  several 
currently  proposed  cutting  planes  for  integer  and/or  concave  programming. 


1 


¥ 

,! 

I 

W.P.  95-71-2  j 


( 

Management  Sciences  Research  Report  No*  279  l 

) 

ON  POLAROID  INTERSECTIONS 

by  ; 
Claude-Alain  Burdet 

i 

April  1972 


-fvjgr^ififTffON  S'f  Affer/tcm1  A 

lor  public  release; 

iii.;tri!',n‘:-on  Unlimited 


This  report  was  prepared  as  part  of  the  activities  of  the  Management  Sciences 
Research  Group,  Carnegie-Mellon  University,  under  Contract  N00014-67-A-0314-0007 
NR  047-048  with  the  U.  S.  Office  of  Naval  Research.  Reproduction  in  whole  or 
in  part  is  permitted  for  any  purpose  of  the  U.  S.  Government. 


Management  Sciences  Research  Group 
Graduate  School  of  industrial  Administration 
Carnegie-Mellon  University 
Pittsburgh,  Pennsylvania  15213 


D  D  C 

EffilcLiilEil'n', 


,  .  Jim  12  1912 

Life?-'  'ITityi 


ABSTRACT 


Polaroid  sets  and  functions  have  been  introduced  as  a  new  tool, 
with  applications  in  non-linear  programming,  particularly  in  quasi-concave 
and  integer  optimization  problems  over  a  linearly  constrained  set  of  feasible 
solutions. 

The  name  polar  programming  applies  to  a  general  class  of  non-linear 
mathematical  programming  problems  which  can  be  solved  by  the  polaroid 
approach.  In  integer  programming  polaroids  yield  non-trivial  extensions  of 
the  intersection  cut  approach. 

This  paper  builds  on  the  properties  of  polaroid  sets  (particularly  com¬ 
plete  convex  polaroids)  and  focuses  on  the  following  intersection  problem: 

Given  a  point  x  belonging  to  the  polaroid  set  P  ,  find 

the  intersection  point  u  of  a  one-dimensional  ray  u  with  the 

•j|(  ^  ^ 

boundary  of  P  :  u  e  (bd  P  ft  u)  . 

We  also  present  a  theoretical  comparison  of  the  relative  merits  of  several 
currently  proposed  cutting  planes  for  integer  and/or  concave  programming. 


On  Polaroid  Intersections 


by  Claude-Alain  Burdet 

0)  Introduction 

Folarold  sets  and  functions  have  been  introduced  in  [ 9] ,  where  we  men¬ 
tioned  areas  of  application  in  non-linear  programming,  particularly  in  quasi¬ 
concave  and  integer  optimization  problems  over  a  linearly  constrained  set  of 
feasible  solutions. 

Polaroids  can  be  used  to  improve  the  efficiency  of  the  algorithms  out¬ 
lined  by  Hoang  Tuy  in  [16]  for  the  maximization  (minimization)  of  a  quasi- 
convex  (concave)  object ive  function  over  a  polyhedral  set.  They  also  improve 
in  a  similar  way  the  modified  version  of  this  approach,  presented  In  f 12]  by 
Glover  and  Klingraan.  Further  algorithms  outlined  in  [8]  can  also  be  extended 
by  using  polaroids  and  in  [10]  we  gave  the  name  polar  programming  to  a  gen¬ 
eral  class  of  non-linear  mathematical  programming  problems  which  can  be 
solved  by  the  polaroid  approach. 

In  integer  programming  polaroids  yield  non-trivial  extensions  of  the 
intersection  cut  approach  of  Balas  [1,  11,  5]  as  indicated  in  [9];  in  [7], 
a  special  type  of  polaroid  is  presented  and  we  show  how  it  brings  under  one 
roof  the  enumerative  approach  of  [5,6]  and  the  extension  method  of  Balas  [2,3]  ; 
it  also  illustrates  a  connection  between  concave  and  integer  programming,  of 
a  nature  different  from  that  of  Rajavaghari  [14] . 


2. 


Consider  Che  following  general  principle  for  solving  the  (arbitrary) 
mathematical  programming  problem: 

maximlxe  g(x),  subject  to  xeP 


where 


x  is  an  n  vector 


P  a  subset  of  3 


g  a  real  valued  function  of  x  e  B 


Method:  Construct  a  (finite)  collection  of  n-dimensional  simplices 
SjCg”  ,  jeJ  such  that 


a)  P  c  S  -  U  S. 


b)  VjeJ  one  has 


xe (S j  0  P) 


g(x)  <  A j 


where  is  a  lower  bound: 

A j  <  g  (x)  ,  for  some  x  e  P 

One  easily  establishes  that  the  optimal  solution  x  of  the 
original  problem  is  delivered  by  (finite)  search:  x  e  [xJ,  jeJ} 
such  that  g  (x)  >  g  (x^)  ,  VjeJ  . 

Particular  examples  of  this  general  method  are  described  in  [8] 
and  [10];  they  can  be  classified  as  follows: 

-  cutting  plane  algorithms  (see,  for  instance,  [16],  [1], 

fill  >  [s],  [6]  ,  [7]) 

-  enveloping  algorithms  (see  [10]) 

"  Partitioning  algorithms  (see  [10],  [16],  [A],  [12]) 

One  of  the  basic  ingredients  in  the  above  general  method  lies  in  the 
construction  of  the  simplex  ,  VjeJ  i.e.,  in  the  determination  of  its 
(n+1)  vertices. 
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this  paper  builds  on  the  properties  of  poiaroid  sets  (particularly  com¬ 
plete  convex  polaroids.  as  defined  in  [9])  which  allow  for  the  construction 
of  simplices  S ^  possessing  the  desired  properties.  We  focus  our  atten¬ 
tion  on  the  determination  of  the  vertices  of  Sj  ,  obtained  here  from  the 
following  intersection  problems 

Given  a  point  x  belonging  to  the  Polaroid  set  P  (for  a 

•k 

definition  see  section  1  below),  find  the  intersection  point  u 

of  a  one -dimensional  ray  u  ,  with  the  boundary  of  p  :  u  e  (bd  P  fl  u)  . 

First  we  outline  the  basLc  methodology  which  leads  to  Polaroid  inter¬ 
section  algorithms.  In  a  second  part  we  sketch  some  typical  examples  of 
Polaroid  functions  and  sets. 

The  interested  reader  is  referred  to  [3,4,7,12,19]  for  application  areas 
and  to  [9]  and  [10]  for  further  aspects  of  the  theory  ct  polaroids. 

1)  Definitions 

In  order  to  make  this  report  self  contained,  we  briefly  reproduce  here 
some  definitions  from  [9]. 

Let  the  Polaroid  function  £  =  £(x,y)  be  real  valued  vith  two  arguments 
x  and  y  which  are  n-vectors;  let  ?  denote  a  closed  set  in  ®.n  . 

•Jf 

Definition  1:  The  Polaroid  set  P  (k)  defined  by  the  Polaroid  function  f 
with  respect  to  the  set  P  and  the  parameter  k  e  f(P,  R°) 

;is 

P*(k)  »  .[y  |  f(x,y)  <  k  ,  VxePl 


4. 


One  has  the 

if 

Theorem  1:  The  Polaroid  set  P  (k)  is  convex  Vk 

iff  the  Polaroid  function  f  *  f<x,y)  Is 
guasi-convex  in  y  ,  for  all  xeP. 

Proof:  Presented  in  [9]  and  omitted  here. 

Although  non-convex  polaroid  sets  P  may  also  be  of  occasional  inter¬ 
est  (6ee  [9]),  we  are  mainly  concerned  here  with  convex  polaroids. 

Definition  2:  Consider  the  point  x  e  P*(k)  and  a  vector  a  e  *n  ;  we 
define  the  ray  u  ,  as  the  following  half-line  originating 
at  x  in  the  direction  (-a)  : 


u  -  {u(X)|u=  x  -  Xa  ,  X  >  0} 


Denote  by  F  *  F  (X)  the  real  vaLued  function 

F_(X)  “max  {f(x,x  -  Xa)} 
xeP 

•ff 

Proposition  1:  Let  X  satisfy 

1)  X*  >  0 

2)  F  (\*)  =  k 


,  *  *  ^ 
then  u  =  u(X  )  lies  on  the  "boundary" of  P  (k)  ,  i.e, 

u  e  bd  P  (k)  =  {y  e  P*(k)  [  f(x,y)  =  k  ,  for  some  xeP} 


if 

Proof:  From  the  hypothesis  2),  3  x  such  that 

^  a  if  if  if 

f(x  ,  x  -  X  a)  =  f(x  ,  u  )  =  k  . 


QF-D 


We  now  briefly  indicate  below  (a  more  detailed  study  can  be  found  in  [9]) 
how  polaroids  find  their  justification  and  application  in  non-linear  optimiza¬ 
tion  problems:  the  following  results  make  it  possible  to  solve  (in  principle) 
a  very  broad  class  of  problems  using  polaroids  and  the  general  method  outlined 
above . 


Definition  3:  The  polaroid  P  (k)  is  called  complete  if  P  c  p  (k)  , 

3b.S.P£«jn„2s  Let  A  “  gtfc)  -  max  g(x),  with  g(x)  -  f(x,x)  .  if  P*<a)  t8  complete 

xeP 

then  every  optimal  solution  x  of  the  problem 

maximize  g(x)  ,  subject  to  xeP 

lies  on  the ’boundary1  of  P  (A)  . 

Prooft  Presented  in  [9]. 

* 

Proposition  2:  Assume  f(x,y)  quasi-convex  in  y  ,  VxeP  ;  then  P  (A)  is 

a  complete  convex  Polaroid  and  any  collection  of  simplices 

S,  ,  jeJ  such  that  S  c  P*(A)  VjeJ  and  P  c  S  =  U  S 

1  j  e  J  j 

yields  a  sufficient  optimality  condition. 

‘k 

Proof:  Convexity  of  P  (A)  was  established  in  Theorem  1  and  completeness 
follows  from  the  choice  of  A  ;  optimality  of  x  then  follows  from 
the  observation  that  f(x,y)  <A  =  g(x)  ,  VxeP  ,  VyeSj  ,  VjeJ 
and  from  the  assumption  PCS.  QED 

The  general  method  described  in  the  introduction  now  becomes  operational 
when  the  simplices  are  constructed  as  subsets  of  the  polaroids  P  (k)  , 

because,  for  k  <  0,  one  has  Sj  c  P*(k)  c  P*(A)  (a  proof  is  given  in  [9]); 
in  practice,  the  value  of  k  is  gradually  increased  (step-wise),  until  the 

optimal  value  A  is  reached.  Thus  the  only  problem  remaining  in  the 
design  of  an  algorithm  is  the  practical  construction  of  the  simplices  . 

2)  The  Intersection  Method 

The  proposition  1  indicates  that  the  determination  of  intersection 
point  reduces  to  the  parametric  linearly  constrained  optimization 
process  described  below: 

Find  X  >  0  such  that  F  (X)  —  max  f(x,x  -  Xa)  <  k  ; 

a  xeP 

if  Pfl(X)  <  k  then  u  =  u(X)  can  be  seen  to  lie  in  the  "interior  of  P*; 

Int  P*(k)  =  [y  |  f (x,y)  <  k  ,  VxeP]  ; 
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i 

I 

i 

i 

i 

j 

I 


if  F#(X)  ■  k  ,  than  u  lias  on  the"boundary" aet  bd  P*(k)  defined  in  sec- 

tion  1  and  it  la  therefore  a  solution  u  to  the  intersection  problem 
Hie  Polaroid  intersection  problem  therefore  can  be  seen  to  contain  the  fol¬ 
lowing  two  parts; 

Pi)  A  Mathematical  programing  part: 

maximize  z(x)  *  f(x,x  -  Xa)  ,  \  given 
subject  to  xcP 

P2)  A  paifiuimt'tHc  problem; 

Find  X*  “  max  {X  I  F  (\)  <  k} 

X>0  a 

In  general  neither  of  these  optimizations  can  be  executed  exactly  and  one  often 
prefers  a  good  approximate  solution  which  is  obtained  in  few  computations. 

Ihis  approximation,  however,  must  be  such  that  the  approximate  intersection 

A.  ^ 

point  u  belongs  to  the  set  P  (k)  . 

For  the  problem  PI)  this  means  that  one  needs  an  upper  bound  solution 
x  such  that 

z(x)  <  It  -  z(x)  =>  F  (X)  <  k  . 
xsP  a 

For  tha  problem  P2)  this  implies  X  <  X  . 

We  now  turn  to  the  important  special  class  of  convex  Polaroid  sets,  there 
the  intersection  problem  plays  a  fundamental  role  in  the  construction 
of  algorithms  (see  [21J). 

Corollary  1.1;  If  tha  polaroid  function  £(x,y)  Is  quasi-convex  in 

y  ,  VxeP  ,  then  VaeRn  the  intersection  function  F  (\) 
i9  quasi-convex  in  \. 


■iHisfcl'miiW  lH  iff— 


7. 


Proof:  By  definition  on*  has  Vac It n  : 

Fa(X)  >  f(*,y  “  x  -  Xa)  ,  VxcP 

Let  y  “  x  -  X  a  ,  y  ”  *  -  X  a  ,  X3=y,X1  +  (l-p,)\2  with  p,c[0,l]  : 

T3  -  x  -  [(iX1  +  (l^)X2]a  -  p.yl  +  (l-n)y2  ; 

for  1*1, 2, 3  one  has  by  hypothesis  Ffl(X1)  >  f(x,yi)  ,  VxeP  ,  with 
Pa(Xi)  “  f^.y1)  ,  for  some  x^P  ; 

hance  Fa(X3)  -  f(x3,y3)  «  f(x3,^y1+  (l~M,)y2)  <  max  {f^.y1),  f(x3,y2)} 

<  max  {f(x1,y1)«  f(x2,y2)}  -  max  {Fa1),  F  a2)} 

w  o  A  J 

Q.E.D. 

Theorem  3:  Let  f(x,y)  be  quasi-convex  in  y  ,  VxeP  ;  and 

let  X*  *»  max  {X  I  F  (X)  <  k} 

X>0  a 

then  Fa(X)  <Fa(X*)  -k  ,  VXc  [0,X*] 

Proof:  By  hypothesis  x  c  P*(k)  ,  implying  f(x,x)  <  k  ,  VxeP  thus 

Fa<0)  <  k  ;  the  quasi-convexity  of  F&  established  in  the  corollary 
1.1  completes  the  proof. 

Q.E.D. 

The  above  theorem  corroborates  a  well-known  property  of  convex  sets 
with  respect  to  their  intersection  by  a  one-dimensional  line  (ray);  the  situation 
Is  complicated  here;  by  the  fact  that  the  function  Fa(\)  is  not  explicitly 
known  but  merely  defined  in  terms  of  the  (arbitrary)  optimization  problem  pi; 
naturally  the  determination  of  an  intersection  point  could  also  be  formulated 
asfone  single  optimization  problem  in  the  variables  x  and  X  ;  however, 
the  independence  of  the  constraints  in  x  and  X  variables  respectively  as 
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well  as  the  separate  properties  of  f  (with  respect  to  x  and  y  respec¬ 
tively)  motivates  a  separate  treatment  of  the  optimization  in  x  and  in  X  ; 
the  result  of  theorem  3  indicates  that  the  parametric  search  problem  P2  can 
be  solved  by  increasing  stepwise  the  ^value  of  X  until  X  (or  X  <  X  ) 

is  reached;  each  step  consistsin  a  (post-)  optimization  of  the  mathematical 

\ 

programming  problem  PI.  The  resulting  algorithm!  becomes  a  non-linear  descent 

algorithm  of  a  particular  type.  \ 

'  »  \ 

Algorithm: 


Step  0:  Let  XQ  =  0  ,  i  =  1  ;  choose  X ^  >  0  such  that 
f(x,x  -  X^a)  <  k  ,  VxeP  ,  i.e.  x  ^  X^  a  »  y^  c 


P*(k) 


Step  1:  Solve  the  (linearly  constrained)  mathematical  programing 


problem 
PI;  max  z 


(J)  —  fi(x,x  -  X  a)  ,  s.t. 

\  I  . 


Step  2: 


2.1s  If  z  =  max  z(x)  =  k  then  jSTOP | 
xeP 


2.2:  lr  z  >  k  ,  then  choose 


Xi+1  e  (X1-1’  V 


2.3;  If  z  <  k  ,  then  choose 


Xi+1  >  Xi 


set  i  -  i  +  1  and  go  to  Step  l. 


>/  8', 


Remarks: 

1)  The  above  algorithmic  principle  merely  represents  an  unpolished  set  of  guide- 

\ 

lines;  depending  on  the  particular  problem  at  hand  it  can  be  refined  to 

i 

increase  computational  efficiency  (see  section  3). 

2)  Clearly  the  above  algorithm  nee (5  not  be  finite;  practical1'’  however,  one 

\  w 

should  remember  that  there1  usually  is  no  need  for  an  accurate  value  of  X  ; 
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i 

thus  I  the  exit  critarion  of  step  2  normally  will  be  made  to  contain  a  (often 
large)  "fuzz"  factor  c  >  0  and  reads: 

I 

2.1:  If  z  e  [k-s  ,  k}  then  STOP  \ 

and  2.3  is  modified  accordingly  to  become 

2.3:  If  x  <  k-s  then  choose  X..,  >  X. 

.  i+1  i 

The  algorithm  then  delivers  an  approximate  intersection  point  u  e  P  (k)  in  a 

finite  number  of  iterations  (this  number  clearly  depends  on  the  quantity  e) . 


3)  Some  Examples  of  Polaroid  Intersections 

\ 

This  section  presents  below  a  list  of  particular  polaroids}  ranked  by 

order  of  imcrdasing  complexity  of  the  corresponding  intersection  problems. 

\ 


3.1)  Bilinear  polaroids:  f(x,y)  «  Ax  +  By  +  y  Cx 


In  this  case  the  polaroid  P  (k)  is  convex  and  the  deter¬ 


mination  of  an  intersection  point 
a  parametric  linear  program. 


*\ 

u  !  is 


obtained  by  solving 

\ 


max  z(x)  =*  B(x-Xa)  +  [A  +  (x-Xa)  C]  x 


\ 


t.t.  xeP 


Problem  PI  is  an  prd inary  L.P.  and  the  increment  of  X  in 
step  2  of  the  algorithm  is  determined  by  the  previous  optimal 
solution  x  of  step  1.  Because  Fft(\)  is  known  to  be  piece- 

•ft 

wise  linear .  the  optimal  value  X  can  be  obtained  by  linear 
algebra  in  a  finite  number  of  Iterations. 

i 

3,1.1. :  For  the  case  where  C  is  a  symmetric  positive  definite 
matrix  Balas  has  first  recognized  in  [3]  the  possibil¬ 
ity  to  use  sets  he  calls  outer-polars  to  generate 

l 

\  \ 

\ 


valid  cutting  planes  In  Integer  programming  with  the 
help  of  parametric  linear  programming  (the  outer- 
polar  is  a  generalized!  polar  set  (see  [15])  and  also 
a  special)  polaroid  set[9]).  The  same  technique  is  used 
by  Bales  and  Burdet  in  [4]  to  solve  quadratic  concave 
programs. 

3.1.2;  For  indefinite  as  well  as  semi -definite  quadratic 

programSj poloroid  cuts  have  been  implemented  success¬ 
fully  in  a  facial  decomposition  schema  for  solving  the 
general  quadratic  programming  problem  [21,22], 

3.U3:  Another  application  for  valid  cuts  in  integer  program¬ 

ming  la  given  in  [7],  there  Polaroid  cuts  are  constructed 
from  a  homogeneous  (A*3®*5!))  bilinear  polaroid  function 
which  is  cantered  at  the  origin  (i.e.at  x)  ;  this 
case  presents  computational  advantages  because  the 
optimization  problem  Pi  and  P2  can  be  combined  in  a 
single  linear  program  (not  parametric);  indeed,  one 
can  verify  by  inspection,  that  the  problem  PI, 

T 

i.e.,  max  -(Xa  Cx)  <  k  , 
xeP 

directly  delivers  the  optimal  value 

*  -1  1  max  (-aT  Cx) 

^  '  ~  k  xeP 

of  problem  P2. 

3.2)  The  previous  bilinear  polaroids  are  a  particular  case  of  the  fol¬ 
lowing  family  (with  parameter  cr): 

f(x,y)  “  g(x)  +  Q'(y-x)1'  F(x)  ,  where  g(x)  =  xT  F (x)  , 

(F  is  a  vector  valued  function). 


As  an  illustration,  consider  the  (quadratic)  case  where 
F(x)  =  Cx  (C  *=  n  by  n  symmetric  matrix)  ; 
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oam  hat  f(x,y)  -  (l-a)  XT  Cx  +  <y  yT  cx 

(For  simplicity,  let  us  assume  C  positive  definite.) 

Atlalf  at  *»  1  :  f(x,y)  “  yT  Cx  ,  bilinear  polaroid  (see  3.1  above) 
hill-  Qf  «  2  ;  f  (x,y)  “  -xT  Cx  +  2yT  Cx  . 

Note  that  ot  F(x)  «  2  Cx  =  Vg(x)  ; 

Thus  we  have  here  an  illustration  (quadratic)  of 
the  general  polaroid  function  f(x,y)  -  g(x)  +  (y-x)T  vg(x)  . 
Bit  Polaroid  set  P*(k)  provides  here  for  an  analytical 
SkftKftctggl of  the  level  set  of  the  maximal  convex 
(concave)  extension  of  the  convex  (concave)  function  g(x) 

ijf-ffa.,ra*£QCt  to  the  set _ P_  ,  which  has  been  characterized, 

in  geometrical  terms,  by  Hoang  TUy  in  [16]. 
at  >  1  :  In  this  case  (which  contains,  in  particular, 
the  previous  one  3.2.2),  the  mathematical  programming  prob- 
le”  P1  18  a  convex  quadratic  programming  problem  m,„r  » 
polyhedral  feasible  set  P  and  it  can  be  solved  by  the 
corresponding  classical  algorithms. 

a  “  ®  !  As  <*  becomes  very  large,  only  points  y  with 

(y-x)  =  0  will  become  tolerable;  the  polaroid 
eet  therefore  simply  becomes  identical  with  the 
original  set  P  . 

3*24-5!  or  <  1  :  Here  problem  PI  is  a  concave  program  (i.e. 

max(min)imization  of  a  convex  (concave)  objective 
function). 

When  0  <  of  <  1  ,  an  illustration  of  the  unifying 
Insight  provided  by  the  polaroid  approach  in  the  fields  of 
concave  and/or  integer  programming  can  be  found  in  [7], 


We  describe  here  a  more  general  situation: 

fSt  ‘l' 

aesume  that  A  “  g(x)  =  max  x  Cx  ;  then,  for  a  small 

x«P 

enough  value  5  >  0  ,  the  polaroid  P*(A)  will  be  the 
half apace (a)  determined  by  the  hyperplane (a)  supporting  P 
at  the  optimal  solutlon(s)  x  .  ihis  is  the  "largest" 

Polaroid  P  in  the  family,  with  the  desirable  property 
tliat  5  ft  Int  P  —  (P  -  bd  P  )  ,  as  shown  in  the  proposi¬ 
tion  3  below,  Larger  values  at  >  at  yield  a  lower  bound  for 
the  intersection  parameter  X  ,  which  can  be  used  as  an 
approximation;  the  intersection  algorithm  could  thus  conceivably  be 
refined  by  changing  the  parameter  a  (starting  with,  say, 
the  L.P.  value  at  ~  1  )  and  post-optimizing  the  (new) 
problem  P2  to  find  a  better  value  for  X  ;  ultimately  this 
process  determines  the  intersection  parameter  X  with  the 

k 

Polaroid  P  (A)  . 

5 

In  integer  programming,  the  outer-domain  theory  (see, 

for  instance, [6]  and  [7]) indicates  that  P*(A)  is  the 

or 

Intersection  of  a  collection  of  halfspaces,  each  belonging 
to  a  feasible  vertex  of  the  unit  cube:  decreasing  the 
parameter  at  (down  from  ot  —  1  )  therefore  corresponds 
to  introducing  additional  integrality  reouirements  into 
the  linear  program  obtained  in  3.1,  see  [7]. 

0=0  :  In  this  case  f(x,y)  -  g(x)  and  is  independent 

k 

of  y  ;  hence  P  is  the  whole  space  whenever 
kef (P,lRn  ), and  P*=0  otherwise.  j 
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*  T  T 

Proposition  3:  Define  Pft(k)  "  (y  J  x  Cx  +  e*(y-x)  Cx  <  k,  VxeP|  ; 

Assume  a  >  0  snd  P  c  leVjcg(x)>  that  is:g(x)  -  xTCx  <  k,  VxeP 

tten  P*  (k)  c  P*  (k)  iff  of.  >  a,  . 

“l  2 

Proof;  Talcs  y«  «  bd  P_  (k)  and  choose  x,cP  such  that 

111  1  *• 

T  x 

f(xJ»y1)  *  \  Cxx  +  oa(yl“xr)wC^l  -  k  ; 

Cx^ 

than  one  has  <y,  -x,  )T  Cx,  >  0  because  0  <  or,  «  • . “  -r  '  » 

1  *  1  <yl-«l>  cx1 


Phare  «**<*,  <  k  |  by  hypothesis. 


Furthermore  y^P^  (k)  iwpll®*  “2  -  “l  *  38  8hown  below: 

xi  Cxi +  «2<yrxi)^Cxi  *  xi  Cxi +  “i(yrxi)Tcxi +  («2-«i)<yrxi)Tcxj  - 

hence  (a2“°j)  is  0  must  hold  true 

since  x^sP  and  x^  Cx^  +  o(1(y1*x1)TCx1  »  k  ,  with  (y^-Xj^^Cx^^  >  0  . 


rP 

On  the  other  hand,  assuming  «2  ~  al  *  8UPPose  3^  ^  > 

then  3  x2  «  P  such  that 


T  T 

x2’Cx2  +  a2(y1-x2)  Cx2  =  >  k; 

_  k'-x,  Cx2 

but  (yi“X.>)  Cx,  >  0  because  0  <  a,  51  - “ — 5 - 

122  1  (yx-x2)  cx2 


Hence  one  obtains  from  y^^  e  PQ^(k)  : 

x2 'Cxj+  «2(y1-x2-)T  Cx2  +  ^yi“x2^T  Cx2  -  k 
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that!*  (or^-o^)  <  0  must  hold  because 

^  0*2  +  a2^yi“x2'  Cx2  >  k  by  hypothesis  and  (yj-x2)T  Cx2  >  0  . 

But  phis  implies  ort  <  «2  which  is  contrary  to  hypothesis;  thus 

y,  «  P*  (k)  .  Q.E.D. 

1  a2 

One  my  conclude  from  proposition  3  that  there  may  be  a  reward  for  solving 
increasingly  difficult  problems  since  for  small  or  the  intersection  prob¬ 
lem  la  of  the  same  order  of  difficulty  as  the  original  (for  instance,  concave 
or  integer)  problem,  one  is  trying  to  solve  with  the  help  of  polaroids. 

Mote  that  parametric  linear  programs  (a  =  1)  however  hold  a  position 
which  seem*  to  make  them  computationally  quite  attractive  in  that  respect. 

This  class  contains  all  the  previous  cases  and  is  charac¬ 
terized  by  functions  f(x,y)  which  are  linear  in  v  .  VxeP  ;  the  con- 
vexity  of  the  corresponding  polaroid  sets  P  (k)  is  readily  verified  (see 
theorem  1). 

Of  course  this  class  may  contain  a  variety  of  different  types  of  optimiza¬ 
tion  problems  PI:  convex,  concave,  discrete,  etc..  ,  .  but  linearity  in 
y  makes  it  easier  to  predict  adequate  increment  for  X  in  the  algorithm. 

One  should  also  note  that  the  convexity  of  P  makes  a  cutting  plane 
approach  possible  for  problems  which  do  not  possess  convex  level  sets 
of  the  objective  function  g(x).  An  example  of  this  type  can  be  found 
in  [21]  and  a  general  discussion  of  polar  programs  is  presented  in  [10]. 


■’•u.  •«*  ,v.'  fl-ti  £S$5tai5*f2l£@j| 
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3.4  A  more  general  class  yet  is  defined  by  polaroid  functions  f(x,y) 

which  are  merely  required  to  satisfy  the  hypothesis  of  theorem  1  in  order 
to  yield  convex  Polaroid  sets.  Here,  except  for  very  special  functions, 
there  are  no  artifices  available  to  render  the  problem  P2  more  easily 
tractable,  and  one  has  to  resort  to  the  stepwise  incremental  method  des¬ 
cribed  in  the  algorithm  to  solve  the  intersection  problem  (depending  on 
the  actual  function  f  it  may  become  more  or  less  difficult  to  determine 
adequate  increments  X.^)  )• 

3.3  One  nay  also  use  polaroid  sets  which  are  not  assumed  in  advance  to  be  con¬ 
vex;  after  solving  the  intersection  problems  Involved  by  a  cutting  plane 
approach,  for  instance,  one  has  then  to  test  the  validity  of  the  cut;  this 
task  is  essentially  different  from  checking  a  quasi-convex  property  of 
£(x,y):  and  in  some  instances,  cuts  may  a  posteriori  be  easily  proved 
valid:  when  generated  from  non-convex  polaroids,  they  may  also  be  deeper 


than  intarsection  cuts  from  convex  outer-domains. 
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